Anomalous Roughening in Experiments of Interfaces in Hele-Shaw Flows 

with Strong Quenched Disorder 
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We report experimental evidences of anomalous kinetic roughening in the stable displacement 
of an oil-air interface in a Hele-Shaw cell with strong quenched disorder. The disorder consists 
on a random modulation of the gap spacing transverse to the growth direction (tracks). We have 
performed experiments varying average interface velocity and gap spacing, and measured the scaling 
exponents. We have obtained f3 ~ 0.50, (3* ~ 0.25, a ~ 1.0, cti oc ~ 0.5, and z ~ 2. When there is 
no fluid injection, the interface is driven solely by capillary forces, and a higher value of /3 around 
P — 0.65 is measured. The presence of multiscaling and the particular morphology of the interfaces, 
characterized by high slopes that follow a Levy distribution, confirms the existence of anomalous 
scaling. From a detailed study of the motion of the oil-air interface we show that the anomaly 
is a consequence of different local velocities over tracks plus the coupling in the motion between 
neighboring tracks. The anomaly disappears at high interface velocities, weak capillary forces, or 
when the disorder is not sufficiently persistent in the growth direction. We have also observed the 
absence of scaling when the disorder is very strong or when a regular modulation of the gap spacing 
is introduced. 

PACS numbers: 47.55.Mh, 68.35.Ct, 05.40.-a 
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I. INTRODUCTION 



The study of the morphology and dynamics of rough 
surfaces and interfaces in disordered media has been a 
subject of much interest in last years Q], and is an active 
field of research with relevance in technological applica- 
tions and material characterization. These surfaces or 
interfaces are usually self-affine, i.e. statistically invari- 
ant under an anisotropic scale transformation. In many 
cases, the interfacial fluctuations follow the dynamic scal- 
ing of Family- Vicsek (F V) M , which allows a description 
of the scaling properties of the interfacial fluctuations in 
terms of a growth exponent (3, a roughness exponent a, 
and a dynamic exponent z. However, many growth mod- 
els [|, |5|, [?], |^], experiments of fracture on brittle 
materials, such as rock |)| or wood [jl0|, [ll], [l2|], as well as 
experiments on sputtering fll3f , molecular beam epitaxy 
PI , and electrodeposition | |15| have shown that the FV 
scaling is limited, and a new ansatz, known as anomalous 
scaling |7, 
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, 17| , must be introduced. The particularity 
of the anomalous scaling is that local and global interfa- 
cial fluctuations scale in a different way, so in general the 
global growth and roughness exponents, (3 and a respec- 
tively, differ from the local ones 0* and a/ oc . 

In experiments of fracture on brittle materials, while 
a is clearly dependent on the material and the fracture 
mode, the measured exponents ol\ oc seem to have a uni- 
versal value, in the range afB. = 0.8 — 0.9 for three- 



*loc 

dimensional fractures, and af® — 0.6 — 0.7 for two 



dimensional fractures. Lopez et 



reported a = 1.2 



in an experiment of fracture in a granite block, while in 
fracture of wood Morel et al. |l(j obtained a — 1.60 and 
a = 1.35 for two different wood specimens. In these ex- 
periments they also obtained (3 — 0.26 and j3* = 0.14. 
An interesting observation from the experiments of frac- 
ture in wood, where tangential and radial fractures were 
studied jl2| , is that the morphology of the interfaces for 
the former were characterized by large slopes, in contrast 
with the smoother appearance of the latter. As pointed 
out by Lopez et al. |9j, the origin of anomalous scaling 
in these experiments seems to be related with these large 
slopes in the morphology of the interfaces. 

In a previous work [ fl9| we presented an experimental 
study of forced fluid invasion in a Hele-Shaw cell with 
quenched disorder. We focused on the limit of weak cap- 
illary forces, and studied different types of disorder con- 
figurations at different interface velocities and gap spac- 
ings. The results obtained showed that the interfacial 
fluctuations scaled with time through a growth exponent 
[3 ~ 0.50, which was almost independent of the exper- 
imental parameters. The values of the roughness ex- 
ponents were characterized by two regimes, a.\ at short 
length scales and «2 at long length scales. For high in- 
terface velocities, we obtained a\ ~ 1.3 and a-i ~ 0, 
independently of the disorder configuration. For moder- 
ate and low velocities, however, different values were ob- 
tained depending on the disorder configuration, and dif- 
ferent scaling regimes were reported. In addition, when 
the disorder was totally persistent in the direction of 
growth, and for sufficiently strong capillary forces, a new 
regime emerged. This new regime must be described in 
the framework of anomalous scaling and was discussed 
in a recent Letter |20 . We obtained the scaling expo- 
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nents for a specific set of experimental parameters, and 
explored the range of validity of the anomaly. We also in- 
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troduced a phenomenological model that reproduces the 
interface dynamics observed experimentally and the scal- 
ing exponents. 

The objective of the present paper is to analyze the 
experimental evidences of anomalous scaling in a sys- 
tematic way, introducing additional techniques for the 
detection of anomalous scaling, such as multiscaling and 
the statistical distribution of slopes. The experimental 
range of anomalous scaling is studied in depth, and we 
investigate the effects of average interface velocity, gap 
spacing, and disorder configuration. We also consider 
the situation in which the interface is solely driven by 
capillary forces, which allows understanding better the 
behavior of the interface at the scale of the disorder and 
its relation with the scaling exponents in particular limits 
of the experiments. 

The outline of the paper is as follows. In Sec. II 
we review the Family- Vicsek and the anomalous scaling 
ansatzs. In Sec. [II we introduce the experimental setup, 
and in Sec. IV we describe the methodology used in data 
analysis. The main experimental results are presented 
and analyzed in Sec. The discussion and final conclu- 



sions are given in Section VI 



II. SCALING ANSATZS 

A. Family- Vicsek scaling 

The statistical properties of a one-dimensional inter- 
face defined by a function h(x,t) (interface height at po- 
sition x and time t) are usually described in terms of the 
fluctuations of h. The global width in a system of lateral 
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where 



size L is defined as w(L,t) = ([h(x,t) — h(t)] 
x is the position, h(t) = (h(x,t)) x , and < > x is an spa- 
tial average in the x direction. In the dynamic scaling 
assumption of Family- Vicsek (FV) @, w(L,t) scales as 



w{L,t)=L a f{L/t l '% 



(1) 



where z is the dynamical exponent and /(it) the scaling 
function 
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const 



if 



if 



u < 1, 
u » 1. 



(2) 



Here a is the roughness exponent and characterizes the 
scaling of the surface at saturation, w(L,t 3> L z ) ~ L a . 
L z characterizes the saturation time, t x ~ L z . It is cus- 
tomary to introduce a growth exponent (3, which char- 
acterizes the growth of the interfacial fluctuations before 
saturation, w(L,t <C L z ) ~ t" . The exponents a, (3, and 
z verify the scaling relation a — z/3. 

An important feature of the FV assumption is that the 
scaling behavior of the surface can also be obtained by 
looking at the local width over a window of size l<L. 
The scaling is then 



w(l,t) 



if 
if 



t -C P, 



(3) 



The roughness exponent a can also be obtained by 
studying the power spectrum of the interfacial fluctu- 
ations, S(k,t) = (H(k,t)H(-k,t)), where H(k,t) = 
XL[' l ( a; !*) — h(t)] exp(ikx) . In the FV assumption, the 
power spectrum scales as 



S(k,t) = k-^+^sikt 1 ^), 
where s is the scaling function 



s(u) 



u i » 
const 



if 
if 



u < 1, 
u > 1. 



(4) 
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B. Anomalous scaling 

As mentioned in the introduction, experimental results 
and several growth models that do not follow the FV dy- 
namic scaling have led to the introduction of a new, more 
general dynamic scaling hypothesis known as anomalous 
scaling jj, |l6|, [Tij| . The origin of this new scaling comes 
from the different behavior of global and local interfacial 
fluctuations. While the global width scales as in the FV 
scaling, the local width follows a different scaling given 
by 



w(l,t) 




if t < I*, 

if P<i« L z , (6) 
if L z < t, 



where (3* = (a — ai oc )/z is the local, anomalous growth 
exponent, and a/ oc the local roughness exponent. It is 
common to define 6 — a — ai oc as the anomalous ex- 
ponent, which measures the degree of anomaly Jl6| , |17[ . 
Therefore, in the case of anomalous scaling, two more 
scaling exponents (3* and ai oc are needed to characterize 
the scaling behaviour of the surface. When 6 = and 
(3* = local and global scales behave in the same way, 
and we recover the usual FV scaling. 

It is common to define a scaling function for the be- 
havior of w(l,t), defining g(u) as gQ/t 1 '*) = w(l,t)/l a . 
From Eq. (o), g(u) is expected to scale as 
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-(a-a loc ) 



if 



if 



U < 1, 
U » 1. 



(7) 



One of the implications of the anomalous scaling is that 
the local width saturates when the system size saturates 
as well, i.e. at times t x , and not at the local time t; ~ 
l z as occurs in the FV scaling. It is also interesting to 
notice that the anomalous scaling can be identified in 
an experiment by plotting the local width w(l,t) as a 
function of time for different window sizes. According to 
Eq. (^) the different curves would show a vertical shift 
proportional to l aioa . 

In terms of the power spectrum, Eq. (H), the anomalous 
scaling follows a new scaling function s[u) given by 



s(u) 



if 



,20 



if 



B«l, 

u > 1. 



(8) 
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plate 




FIG. 1: Sketch of the experimental setup, a) Side view, b) 
Top view, c) Schematic representation of the oil-air interface 
and the disorder pattern. 



For times t^> t x the power spectrum scales as S(k, t) ~ 

k -(2a loc + l) t 29/z and nQt simply ag fc-(2a+l) of the ord j_ 

nary scaling. 



In general, it is usual to differentiate between two 
types of anomalous scaling: super-roughness and intrin- 
sic anomalous scaling. Super-roughness is characterized 
by aioc = 1 and a > 1, although the structure factor fol- 
lows the ordinary scaling given byEqs. §) and (|). Thus, 
for super-roughness, the power spectrum provides infor- 
mation about the global roughness exponent a. The in- 
trinsic anomalous scaling follows strictly the scaling given 
by the scaling function of Eq. (||) . This implies that the 
decay of the power spectrum with k gives only a mea- 
sure of cti OCl not of a. Moreover, the particular temporal 
dependence of the scaling of the power spectrum allows 
identifying the existence of intrinsic anomalous scaling by 
the presence of a vertical shift of the spectral curves at 
different times. One has to be cautious with this partic- 
ularity, because the observation of a vertical shift, how- 
ever, does not imply necessarily the existence of intrinsic 
anomalous scaling. 



The different forms of anomalous scaling are included 
in a generic dynamic scaling ansatz recently introduced 
by Ramasco et al. 21 . 



III. EXPERIMENTAL SETUP AND 
PROCEDURE 

The experimental setup used here has been described 
previously in Refs. M. The setup (Fig. |l|) consists 
on a horizontal Hclc Shaw cell which contains a pattern 
of copper tracks placed over a fiberglass substrate. The 
tracks are d=0.06 ± 0.01 mm high and are randomly dis- 
tributed along the lateral direction x without overlap, 
and filling 35 % of the substrate. The lateral size of the 
unit track is 1.50 ± 0.04 mm, but wider tracks are ob- 
tained during design when two or more unit tracks are 
placed one adjacent to the other. The separation be- 
tween the fiberglass substrate and the top glass defines 
the gap spacing b, which has been varied in the range 
0.16 < b < 0.75 mm. A silicone oil (Rhodorsil 47 V) 
with kinematic viscosity v = 50 mm 2 /s, density p — 998 
kg/m 3 , and surface tension oil-air a = 20.7 mN/m at 
room temperature, is injected at one side of the cell, and 
displaces the air initially present. The evolution of the 
oil-air interface is monitored using 2 CCD cameras. A 
typical example of the morphology of an interface and its 
temporal evolution is shown in Fig. 0. 

The oil is injected at constant volumetric injection rate 
Q by means of a syringe pump. The presence of the dis- 
order, however, makes the cell inhomogeneous in the z 
direction, so that Q is not exactly proportional to the av- 
erage interface velocity v measured on the images. Nev- 
ertheless, since the free gap b — d is always comparable 
to the total gap 6, the two-dimensional interface velocity 
v is still well characterized. On the other hand, to make 
the comparison between different experiments easy, we 
will measure v in terms of a reference velocity V = 0.04 
mm/s. 

Although the disorder made of tracks (which will be 
called T 1.50, the numeric value indicating the width of 
the unit track in mm) has been the main disorder pattern 
used in the experiments, other disorder configurations 
have been used. These configurations are characterized 
by a random distribution of squares with a lateral size 
of 1.50 or 0.40 mm, and will be called SQ 1.50 and SQ 
0.40 respectively. A full description and characterization 
of these disorders, plus more detailed information about 
the experimental setup and procedure can be found in 
Ref. M. 



IV. DATA TREATMENT 

As explained in Sec. II, we study the roughening pro- 
cess of the interfaces through the interfacial rms width 
w(l,t) and the power spectrum S(k,t). The profile 
h(x,t), obtained from the analysis of the images, con- 
tains 515 points corresponding to a lateral size of 190 
mm. Because the oil tends to advance at the side walls 
of the cell, we have a distorsion of the front that pene- 
trates a distance p from each side wall. This penetration 
length depends on average interface velocity and is spe- 
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(a) i 



FIG. 2: Close-up views of the oil-air interface with disorder T 
1.50. The first three pictures are early stages, while the last 
one corresponds to saturation. The position of the copper 
tracks is shown in the bottom box. The dotted line shows the 
position and shape of the interface profile at t = s. The 
experimental parameters are v — 0.08 mm/s, b — 0.36 mm. 



cially sensitive to variations in gap spacing. An example 
is presented in Fig. ||(a), which shows the effect of the 
distortion on an interface with experimental parameters 
v = 0.08 mm/s, b = 0.75 mm and disorder SQ 1.50. To 
evaluate the effect of the distortion, we have measured 
the penetration length p at different gap spacings b for 
two disorder configurations, SQ 1.50 and T 1.50. The 
results are presented in Fig. 0(b), and show that p in- 
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FIG. 3: (a) Example of an interface with strong side wall 
distortions. The dashed lines represent the regions cut at the 
end of the interface, and the solid line the interface finally 
studied. The experimental parameters are v = 0.08 mm/s, 
b — 0.75 mm and disorder SQ 1.50. (b) Penetration length p 
as a function of the gap spacing b for SQ 1.50 (dots) and T 
1.50 (triangles). 



creases quickly with gap spacing. For b = 0.75 mm, the 
largest value studied, p is about 8 mm. Notice also that 
the presence of a columnar disorder reduces significantly 
the effect of the distortion for the smallest gap spacings. 
To minimize the effect of the distortion we have cut 8 mm 
from both ends of all interfaces, as shown in Fig. ||(a), 
reducing them to a final lateral size of 174 mm (470 pix- 
els). For data analysis convenience (particularly in the 
computation of the power spectrum), we have extended 
the number of points to 512 using linear interpolation. 
Finally, we have forced periodic boundary conditions by 
subtracting the straight line that connects the two ends 
of the interface. This procedure, well documented in the 
literature of kinetic roughening p2] j, eliminates an arti- 
ficial overall slope -2 in the power spectrum due to the 
discontinuity at the two ends p3| . 

Because the initial stages of roughening are so fast, the 
interface profile at t = s is not a flat interface, but a pro- 
file with an initial roughness u>(0), as can be observed in 
Fig. |. This roughness, on one hand, does not permit to 
observe the growth of the small scales at short times; on 
the other, w(0) is different from one disorder configura- 
tion to another, which makes the early stages strongly de- 
pendent on the particular realization of the disorder. To 
eliminate these effects, we will always plot the subtracted 
width W[l, t), defined as W(l, t) = (w 2 (l, t)-w 2 (l, 0)) 1/2 , 
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where the brackets indicate average over disorder real- 
izations. This correction is an standard procedure in the 
field of growing interfaces, and is reported in both exper- 
iments and numerical simulations Jl], [u| |24|] . 

Another aspect of interest concerning the analysis of 
our experimental data is the various methods available 
in the literature to obtain the growth and roughness ex- 
ponents. In this work we have obtained the growth expo- 
nents (3 and (3* through the scaling of W(l, t) with time. 
We have also used an alternative way to measure 0*, 
proposed by Lopez ]l8| , in which the growth of the mean 
local slopes, p(t) = ((V/i) 2 ) 1 / 2 , is analyzed. Usually p(t) 
scales with time as p(t) ~ i K , with n < for the ordinary 
scaling and k > for anomalous scaling. In this last 
case, k is identified with the local growth exponent (3* . 
The local roughness exponents oti oc have been obtained 
through the analysis of the power spectrum S(k,t). The 
results obtained have been checked using other methods, 
such as the dependency of the height-height correlation 
function C(l, t) or W(l, t) with the window size I, in both 
cases scaling as l a '° c . We have obtained, within error 
bars, similar results to those from S(k,t). Finally, in the 
framework of intrinsic anomalous scaling, the only direct 
way to measure the global roughness exponents a, and 
in addition z = a/f3, is through the analysis of C(L,i) 
or W(L, t) with the system size L. Both quantities scale 
as L a . Because it is not practical to perform experi- 
ments varying the system size, we have not measured 
global roughness exponents directly. Instead, they have 
been obtained from the collapse of S(k, t) or W(l, t). No- 
tice that information about the value of a is contained 
in both the vertical displacement of the spectral curves 
with time (Eqs. (Q) and (||)) and the vertical shift of the 
W(l,t) curves for different window sizes (Eq. (||)). The 
error bars that we give in the final numeric value of an 
specific scaling exponent takes into account the various 
methods of analysis and the error in the different col- 
lapses. 
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EXPERIMENTAL RESULTS AND ANALYSIS 



A. Evidence of anomalous scaling 

The first evidence of anomalous scaling in these ex- 
periments was presented in our previous work |20|] and 
corresponded to the experimental parameters v = 0.08 
mm/s (2V) and b = 0.36 mm. Here, for completeness, we 
will show similar experimental results for a different set 
of parameters, namely v = 0.04 mm/s (V) and b = 0.36 
mm. 



1. Interfacial width 

Fig. | (a) shows the results obtained for W{1, t) for pro- 
gressively smaller window sizes. We obtain clear expo- 
nents (3 — 0.52 ± 0.04 (corresponding to I = L) and 



FIG. 4: Experimental results for b — 0.36 mm, v = 0.04 
mm/s (V) and disorder T 1.50. a) Analysis of W(l,t) for 
different window sizes, b) Analysis of the mean local slopes 
as a function of time, c) Temporal evolution of the power 
spectrum. The vertical line in c) indicates the spatial scale of 
the disorder. 



13* = 0.35 ± 0.05 (for I = L/256). The fact that (3* is 
clearly different from zero is a robust evidence of anoma- 
lous scaling. Moreover, the behavior of the different 
curves is in good agreement with the anomalous scal- 
ing ansatz of Eq. (||), i.e. the plots for different window 
sizes are vertically shifted. Notice also that the satura- 
tion takes place at t x (the saturation time of the full 
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system) for all window sizes. 



2. Local slopes 

To check the reliability of the exponent (3* we have also 
studied the scaling of the mean local slopes with time. 
The result is presented in Fig. [|(b), obtaining j3* = 0.38± 
0.05, in agreement with the result of (3* using W(l = 
L/256,i). 



3. Power spectrum 

The analysis of the interfacial fluctuations through the 
power spectrum is presented in Fig. ||(c). The first im- 
portant observation is that the spectral curves are shifted 
vertically, the shift being maximum at short times, reduc- 
ing progressively as time increases, and disappearing at 
saturation (t = t x ~ 2000 s). The presence of this ver- 
tical displacement of the spectral curves is a signature 
of anomalous scaling. The evolution of the power spec- 
trum gives two clear power law dependencies, at short 
times with slope —3.4 ±0.1, and at saturation with slope 
— 2.1 ± 0.1. The former is a transient regime whose ori- 
gin is the strong dominance of the capillary forces at very 
short times. The latter provides the local roughness ex- 
ponent, ai oc = 0.55 ± 0.10 in this case. 

To get the values of a and z we have collapsed the plots 
of Fig. ^(c) in accordance with the scaling function of Eq. 
(§), i.e. plotting S(k, t)k 2a+1 as a function of kt x / z . This 
collapse, shown in Fig. ^|(a), refines the numerical values 
of the experimental results and gives the values of a and 
z. We obtain in this way: 

(3 = 0.50 ± 0.04, (3* = 0.30 ± 0.04, 

a = 1.1 ± 0.1, a ioc = 0.5 ± 0.1, z = 2.2 ± 0.2. (9) 
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FIG. 5: a) Collapse of the power spectrum of Fig. ^(c) for 
t > 10 s using the scaling function of Eq. (|^). b) Collapse of 
the W(l, t) plots of Fig. 0(a) , using the scaling function of Eq. 
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An alternative collapse that confirms the validity of 
these results can be performed using the scaling function 
defined in Eq. ([?]). This collapse is shown in Fig. |B|(b). 



B. Experimental range of anomalous scaling 



4- Multiscahng 

The presence of anomalous scaling can also be iden- 
tified by checking for multiscahng Q. For this pur- 
pose, we have studied the scaling of the general- 
ized correlations of order q, of the form C q (l) = 
{(h(x + l,t) -h{x,t)Y} 1 /i - l a ", where a q are local 
roughness exponents. The results obtained are shown 
in Fig. ||, which gives the values of a q at different or- 
ders q. We get a progressive decrease of the measured 
exponent, varying from oli ~ 0.55 to a§ ~ 0.22, which 
confirms the existence of multiscahng in the experiment 
and, therefore, of anomalous scaling. 



Once the existence of anomalous scaling is well charac- 
terized for a specific set of experimental parameters, we 
are now interested in analyzing the range of validity of 
the anomaly. We have performed a series of experiments 
varying v and fixing the gap spacing to b = 0.36 mm and, 
on the other hand, a series of experiments varying b and 
fixing the average interface velocity to v = 0.08 mm/s 
(2V). Because we can identify the existence of anoma- 
lous scaling by looking at the behavior of the smallest 
scales, i.e. the existence of a (3* ^ 0, we present the 
results by plotting the local growth exponent (3* for the 
different experiments. 
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FIG. 6: Multiscaling analysis through the computation of 
the q-order height-height correlation function for interfaces 
at saturation. The experimental parameters are v — 0.04 
mm/s, b = 0.36 mm and disorder T 1.50. 



FIG. 7: W(t) curves for a window size / = L/256 and different 
interface velocities. The gap spacing is b — 0.36 mm in all 
cases, and V = 0.04 mm/s. The curves are vertically shifted 
for clarity. 



1. Interface velocity 

The results of (3* at different v are presented in Fig. |?|. 
In a wide range of interface velocities, from V to 10 V, 
the measured exponent (3* is robust and not very sensi- 
tive to changes in the velocity, varying from 0.35 ± 0.04 
for velocity V to 0.28 ± 0.04 for velocity 10 V. For ve- 
locities v > 15V, however, the measured value of (3* 
quickly decays to a value close to zero, indicating that 
over a threshold velocity (about 13V for b = 0.36 mm) 
the anomaly disappears. 

In terms of the power spectrum, when the velocity is 
varied from V to 15V, we always observe a vertical shift 
of the spectral curves, and the power spectrum at dif- 
ferent times can be collapsed to get the set of scaling 
exponents. Above the previous threshold velocity, the 
power spectrum follows the ordinary scaling, i.e. with- 
out vertical displacement. 



2. Gap spacing 

The results of /?* at different b are presented in Fig. 
||. It is clear from the different plots that the exponent 
/?* is rather sensitive to variations of the gap spacing. 
Three different behaviors can be observed: i) For small 
values of b (around b = 0.16 mm), where capillary forces 
are dominant, we get (3* = 0.52 ± 0.04, which is identi- 
cal to the corresp ond ing value of (3. This similitude is 
ii) For moderate values of b (in 



discussed in Sec. VE 



the range 0.23 < b < 0.56 mm) we get values of /3* that 
vary from 0.35 ± 0.05 for b = 0.23 mm to 0.25 ± 0.04 for 
b = 0.56 mm. This range of gap spacings is character- 



ized by strong (but not dominant) capillary forces, and 
is the region where the anomalous scaling can be fully 
characterized, iii) For large values of b (b > 0.60 mm) 
we get (3* = 0.10 ± 0.04 for b = 0.75 mm, and (3* is 
nearly zero for larger gap spacings. In this regime capil- 
lary forces are not sufficiently strong, and the anomalous 
scaling disappears. The analysis of (3* through the scal- 
ing of the mean local slopes has given identical results 
within error bars, although the sensitivity of this analy- 
sis to short-time transients makes it less reliable. 

Concerning the power spectrum, in the range 0.23 < 
b < 0.56, the vertical shift of the spectral curves is 
present, and a regime characterized by anomalous scal- 
ing can be described. For b = 0.16 mm, although the 
spectral curves are vertically shifted, an scaling regime 
is not identifiable at all. The power spectra can be col- 
lapsed, but not reliable a, a; oc , and z exponents can be 
obtained. For b > 0.6 mm, a vertical shift in the power 
spectra is not observed, and the scaling of the interfacial 
fluctuations is best described in the framework of ordi- 
nary scaling. 



3. Disorder properties 

Other modifications that we have carried out to study 
the range of validity of the anomalous scaling are the 
track width and the persistence of the disorder along the 
growth direction. Fig. ^ shows the variation of the local 
interfacial fluctuations for four different disorder config- 
urations. The gap spacing and velocity have been fixed 
to b = 0.36 mm and v = 2V respectively. For T 0.40 
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FIG. 8: W(t) curves for a window size I — L/256 and different 
gap spacings. The interface velocity is v = 2V = 0.08 mm/s 
in all cases. The curves are vertically shifted for clarity. 



FIG. 9: W(t) curves for a window size I = L/256 and different 
disorder configurations with progressively shorter persistence 
in the growth direction, b = 0.36 mm and v — 2V = 0.08 
mm/s in all cases. The curves are vertically shifted for clarity. 



we get a similar result from that obtained with T 1.50, 
measuring (3* — 0.29 ± 0.05. However, we have evidences 
that the anomaly decreases as we reduce the width of the 
basic track, and it disappears when the width of the ba- 
sic track is < 0.20 mm, due to the fact that the in-plane 
curvature over the copper track becomes similar to the 
curvature in the z direction. 

Next, we study different disorder configurations with 
progressively smaller persistence in the growth direction. 
For SQ 1.50, although the saturation time for the full sys- 
tem is at t x — 1000 s, we observe that the local interfacial 
fluctuations saturate at much shorter time, which is more 
characteristic of the ordinary scaling. However, we still 
observe traces of anomalous scaling at early times, for 
which (3* = 0.26 ± 0.04. This is somehow not surprising 
taking into account that, for this disorder, the first array 
of squares is equivalent to T 1.50, but only 1.50 mm long. 
Thus, at very short times, the growth follows the same 
behavior as tracks. With a driving velocity 0.08 mm/s, 
the interface takes about 10 — 20 s in crossing the first 
array of squares, which is the same time for which we 
observe the presence of anomalous scaling. For SQ 0.40, 
with t x ~ 300 s, we measure (3* = 0.05 ± 0.04, which is 
indicative of the absence of anomalous scaling. 



4- Characterization by the statistical distributions of slopes 

We have also considered an alternative analysis to con- 
firm the results of anomalous scaling. In a recent work, 
Asikainen et al. |2(| showed that the presence of anoma- 
lous scaling can be associated to the presence of multi- 
scaling on one hand, and to high slopes in the morphology 



of the interfaces on the other. According to these au- 
thors, the distribution function of these slopes, P(Ah), 
with Ah = \h(xi+i — h x \, scales following a Levy dis- 
tribution, P(Ah) ~ A/i -7 , where < 7 < 2; while for 
ordinary scaling P(Ah) follows a Gaussian distribution. 
We have checked this prediction with our experimental 
results. Fig. 10(a) compares the shape of three different 
interfaces with different degrees of anomalous scaling. In 
all cases v — 2V and b — 0.36 mm. All the interfaces are 
in the saturation regime. For T 1.50 (the case most rep- 
resentative of anomalous scaling) the interface presents 
high slopes at the edges of the copper tracks. For T 0.40 
some regions present high slopes, but in general the inter- 
face looks smoother. This case shows anomalous scaling, 
but the degree of the anomaly is much weaker than the 
previous case. For SQ 1.50 the interface is very smooth 
and the interfacial fluctuations are best described with 
the FV scaling, with only traces of anomalous scaling at 
very short times. The interfaces in these three cases are 
similar to those reported in fracture of wood by Eng0y 
et al. in Fig. 1 of Ref. jL2|. In particular, the interfaces 
for T 1.50 or T 0.40 look like interfaces obtained from 
tangential fracture in wood, both characterized by high 
slopes and anomalous scaling. In the same line, the inter- 
faces for SQ 1.50 are similar to those obtained from radial 
fracture. These similitudes indicate that the role of the 
wood fibers is probably equivalent to the role of the cop- 
per tracks in our system as the origin of the anomalous 
scaling. 

The distribution functions for the different interface 
morphologies are shown in Fig. |l^(b) and (c). For T 
1.50 we have two different regimes. Up to Ah ~ 2 mm 
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(dashed line in the plot) we have the contribution of the 
slopes corresponding to the structure inside tracks. For 
Ah > 2 mm (solid line), the contribution to the distri- 
bution function arises from the edges of the tracks. This 
last contribution is the one responsible for the anomalous 
scaling, and scales with Ah with an exponent —0.6 ±0.1. 
A similar result is observed for T 0.40. In this case, if we 
do not consider the contribution of the structure inside 
tracks (up to 1 mm), we get a power law dependency 
with exponent —2.0 ± 0.2. The fact that the anomaly is 
weaker for T 0.40 compared to T 1.50 is illustrated in this 
case by the larger value of the slope. If we compare the 
distribution functions of T 1.50 and SQ 1.50 we observe 
that, in the latter case, there is contribution only from 
small slopes. P(Ah) would scale in this case with an ex- 
ponent —3.6, far from the values expect for anomalous 
scaling. 



C. Experiments at Q=0 

In order to understand the origin of the anomalous 
scaling, the previous results indicate that a more detailed 
investigation about the role of the capillary forces at the 
scale of the disorder has to be carried out. For this rea- 
son we have performed experiments at Q = 0, where 
intcrfacial growth is driven solely by capillary forces. 

The experiments at Q = have been performed by 
initially driving the interface at high velocity until it has 
reached the center of the cell. Then the pump is stopped 
and the experiment started. We have investigated five 
gap spacings, in the range 0.16 <b< 0.75 mm, and char- 
acterized the interface velocity over copper tracks and 
fiberglass channels. For sufficiently small gap spacing, up 
to b ~ 0.36 mm, the morphology of the interface adopts 
a finger-like configuration, which ultimately causes that 
the interface pinches-off at long times. For this reason 
we have reduced the duration of the experiments to 500 
s for all gap spacings. 

Fig. |ll] shows a comparison of the interfaces at the 
end of the experiments with different gap spacings. The 
horizontal line represents the average position of the in- 
terface at t — 0. For b — 0.16 mm the points of the inter- 
face over the different tracks move almost independently, 
with weak coupling between neighboring tracks: the oil 
over copper tracks always advances, while the oil over 
fiberglass channels, except for some particular locations, 
always recedes, driven by the different capillary forces as- 
sociated with the different curvature of the meniscii over 
copper tracks and fiberglass substrate. As we increase 
the gap spacing, the coupling in the motion of the inter- 
face over neighboring tracks is progressively more intense. 
For 6 = 0.36 mm the advance over a specific copper track 
is coupled with its neighbors. The oil advances or recedes 
in a given track or channel depending on the particular 
realization of the disorder. Only the widest copper tracks 
or fiberglass channels have an independent motion. For 
b = 0.56 mm the coupling reaches larger regions of the 




20 40 60 80 100 120 140 160 180 

x (mm) 
(b) (c) 




1 10 1 10 

Ah (mm) Ah (mm) 



FIG. 10: (a) Samples of interfaces for different disorder con- 
figurations and identical experimental parameters v — 0.08 
mm/s and b — 0.36 mm. (b) Distribution function P(Ah), 
with Ah = \h(xi+i — h x \, for local slopes, comparing T 1.50 
and SQ 1.50. (c) The same for T 0.40. See text for details of 
the power law fits. 



cell, and for b — 0.75 the capillary forces are not longer 
sufficient to allow a large deformation of the interface. 

It is also noticeable that the average position of the 
interface does not keep at rest. Indeed, it advances in 
the same direction as the oil over copper tracks. This is 
a consequence of the three-dimensional nature of the ex- 
periment. The total mass of fluid in the cell is conserved, 
but not the area measured on the images. However, the 
measured velocity of advance of the interface is very low. 
For b = 0.36 mm, for example, the velocity is about V/10. 

The first analysis that we have performed with the 
experiments at Q = is studying the growth of interfacial 
fluctuations with time. Fig. [12| shows the W(t) curves for 
four different gap spacings. We observe power laws with a 
growth exponent (3 that decreases from 0.65±0.04 for b = 
0.16 mm to 0.50 ± 0.04 for b = 0.75 mm. Clearly, when 
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FIG. 11: Example of interfaces at Q = with different gap 
spacings and disorder T 1.50 mm. All the interfaces corre- 
spond to t = 500 s. The horizontal dotted line indicates the 
position of the interface at t — s. 



the capillary forces are sufficiently strong to allow for 
both the advance of the oil over tracks and the recession 
of the oil over fiberglass channels, we get a larger value 
of the growth exponent than in all other cases. The same 
behavior is observed for j3* , varying from 0.65 ± 0.04 for 
b = 0.16 mm to 0.15 ± 0.05 for b = 0.75 mm. 

Concerning the behavior of the power spectrum, the 
fact that the experiments do not reach a saturation 
regime makes the analysis of the power spectrum incom- 
plete. Qualitatively, we have observed a vertical shift of 
the spectral curves as time increases for gap spacings up 
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FIG. 12: Log-log plot of W(t) for experiments at Q — with 
four different gap spacings and disorder T 1.50 mm. The 
straight line with slope 0.65 corresponds to a fit of the ex- 
periments with b — 0.16 mm and 0.36 mm. The straight line 
with slope 0.50 is a fit to the experiment with b = 0.75 mm. 
Each curve is an average over two disorder realizations, with 
one run per disorder realization. 



to b = 0.56 mm. For b = 0.75 mm the vertical shift is 
not present, confirming the absence of anomalous scaling 
for values of b > 0.60 mm for any interface velocity. 

The second analysis that we have performed has been 
a detailed study of the motion of the interface over cop- 
per tracks or fiberglass channels separately. In Fig. [l3| 
we represent the variation of the mean height < h > for 
points of the interface that are advancing over copper 
tracks. For gap spacings in the range 0.16 < b < 0.36 
mm, the mean height is a power law of time with an 
exponent 0.50 ± 0.03. For larger gap spacings, the expo- 
nent is progressively smaller, and reaches 0.32 ± 0.03 for 
b = 0.75 mm. We have hints that the exponent tends to 
zero as the capillary forces become even weaker. A sim- 
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FIG. 13: Log-log plot of the average height as a function of 
time for the interface points advancing over copper tracks, 
for different gap spacings. Q = in all cases. The curves are 
vertically shifted for clarity. 



ilar analysis can be carried out for the fiberglass chan- 
nels. In this case we get the same behavior in the range 
0.16 < b < 0.36 mm, with an exponent 0.50 ± 0.03. For 
b > 0.60 mm not clear scaling can be identified. Finally, 
we have also studied the motion on individual tracks, in 
order to know the role of track width. Considering tracks 
whose motion is not coupled with its neighbors we have 
obtained that tracks of different width show the same be- 
havior < h >^ t - 5 , but the different curves are displaced 
a factor proportional to the track width. 

From these studies we conclude that it is possible to 
describe the motion over a single copper track, v\_, or 
fiberglass channel, v % _, by the simple law 

v i ±{t)=A i ± (b,D i ±)t- 1 / 2 for t^ls, (10) 

where A l ± (b, D l ± ) is a prefactor that depends on the gap 
spacing b and the individual track or channel width D % . 
Averaging out the prefactor A l ± over different track or 
channel widths we obtain a characteristic maximum ve- 
locity vm± for a given gap spacing. Eq. ([[o|) becomes 
now 

v±{t) = v M ±(b) r 1/2 for t^ls. (11) 

In general vm + > vm- due to the different gap spacing 
on copper tracks and on fiberglass channels. This differ- 
ence, however, is progressively smaller as b increases. For 
b > 0.36 mm it is reasonable to consider vm + — «m_ = 
vm- The time t = 1 s, determined from systematic ex- 
perimental observations, is a typical time in which the 
velocity jumps from the average interface velocity out of 
the disorder to Vm inside the disorder. 



D. From Q = to low injection rates 

Eq. ( pr| ) gives the velocity over copper tracks or fiber- 
glass channels when Q — (no injection). In the presence 
of injection, when the interface is driven at an average 
interface velocity n / 0, the behavior observed experi- 
mentally is modified in the form 

v±{t)=v±{vM-v)t- 1 ' 2 , (12) 

where % = and v+ and V- are the velocities 

averaged over copper tracks and fiberglass channels re- 
spectively. Thus, when the interface gets in contact with 
the disorder, the velocity on copper tracks jumps to a 
maximum vm in a characteristic time of 1 s and then 
relaxes as £ -1 / 2 to the nominal velocity. On fiberglass 
channels, the velocity is zero or even negative at short 
times, and as time goes on, the velocity increases until 
it reaches the nominal velocity at long times. Fig. || il- 
lustrates this behavior. At short times, t ~ 10 s, all the 
oil over copper tracks is advancing, while the oil on the 
widest fiberglass channels is receding. At t ~ 50 s the 
oil on fiberglass channels has stopped receding. At later 
times, the velocity on different points of the interface 
tends asymptotically to the nominal velocity, reached at 
saturation. 

At low injection rates, the external velocity v can be 
tuned appropriately to have an initial regime in which the 
oil advances over copper tracks and recedes over fiber- 
glass channels, followed by a regime in which the oil ad- 
vances at all points of the interface. This leads to two 
scaling regimes for W(l,t), at short times with a growth 
exponent ~ 0.65, which corresponds to a regime dom- 
inated by strong capillary forces, and at longer times 
with a growth exponent ~ 0.50, which corresponds to a 
regime where capillary and viscous forces are better bal- 
anced. We have illustrated this behavior by performing 
experiments at low v. In particular, we have studied ve- 
locities V/2 and V/10 using a gap spacing b — 0.36 mm 
and disorder T 1.50. The results are presented in Fig. 
|l4| . The difficulty of these experiments is that the inter- 
face pinches-off before it reaches saturation, so we can 
investigate only the growth regime. For V/2 we observe 
the two power law regimes described before. For V/10, 
only the first regime can be observed before the inter- 
face breaks. These results are qualitatively in agreement 
with the analysis of Hernandez-Machado et al. [^5| , who 
predicted two scaling regimes, (5\ — 5/6 at short times 
(surface tension in the plane dominant) and /?2 = f/2 
at long times (viscous pressure dominant). However, the 
presence of anomalous scaling and the complicated role 
of capillary forces in our experiments makes it difficult 
a direct comparison with the predictions of Hernandez- 
Machado et al.. 

Another result than can be obtained from Eq. ( fi"2| ) is a 
prediction for the saturation time t x as a function of the 
average interface velocity v. We can estimate t x for a 
given v and gap spacing b by considering that saturation 
takes place when all tracks reach the nominal velocity 




FIG. 14: log-log plot of W(t) for experiments at low velocity 
with disorder T 1.50 mm. The end of the experiments is the 
moment at which the interfaces pinch-off. The gap spacing 
is b — 0.36 mm. Each curve is an average over two disorder 
configurations, with one run per disorder configuration. 



FIG. 15: Prediction of the dependency of t x with v. Main 
plot: for a given average velocity v, the system saturates when 
the average velocity on copper tracks decays to v. Inset: a 
log- log plot of the t x obtained by the previous method (filled 
dots) as a function of v gives a straight line with slope —2. The 
triangles are the experimental results for T 1.50 and b = 0.36 



v. This is not strictly exact due to the coupling between 
neighboring tracks, but it gives a good estimate of the 
dependency t x (v), as shown in Fig. 15, For a given ex- 



periment with velocity Vi < vm the saturation takes place 
when the function vm t^ 1 / 2 crosses the nominal velocity. 
Thus, as shown in the inset of the figure, t x is expected 
to scale as t x ~ v~ 2 . This result is confirmed by our ex- 
periments with disorder T 1.50 and b — 0.36 mm (open 
triangles in the inset). 



E. Experiments with a regular modulation of the 
gap spacing 

In order to analyze the importance of the coupling 
between neighboring tracks, we have performed experi- 
ments with a regular modulation of the gap spacing, con- 
sisting on copper tracks 6 mm wide separated by fiber- 
glass channels 6 mm wide. Interface velocity and gap 
spacing are set to v = 0.08 mm/s and b — 0.36 mm 
respectively. The regular modulation and the temporal 
evolution of the resulting interfaces are presented in Fig. 
|l6|(a). Due to the symmetry of the regular pattern there 
is no growing correlation length along the direction of the 
interface. For this reason, a dynamic scaling of the inter- 
facial fluctuations cannot even be defined. The growth of 
the interfacial width is due to the independent, equiva- 
lent growth of interface fingers over copper tracks, which 
follow the t^ 1 ! 2 behavior (Eq. Q). The system satu- 
rates when the fingers reach the nominal velocity. This is 
shown in Fig. [l6|(b), which shows the evolution of W(l, t) 
for the global \T— L) and local (I = L/256) scales. The 
figure reveals that both scales grow in the same way. The 



absence of dynamic scaling is also confirmed by the power 
spectrum, shown in Fig. |l6|(c). We obtain a nearly flat 
spectrum, with a large peak at the value of k that corre- 
sponds to the spatial periodicity of the underlying regular 
pattern. We also observe a vertical shift of the spectra 
at different times, caused by the progressive elongation 
of the oil fingers over copper tracks during growth. 

The previous experiments with a regular modulation 
of the gap spacing illustrate the need of a growing cor- 
relation length along the interface to reach a dynamic 
scaling regime. Similarly to these experiments, the lat- 
eral growth of the correlation length can also be inhibited 
when the modulation of the gap spacing is random, pro- 
vided that the experimental parameters are such that the 
motion over tracks is not coupled and, on each track, fol- 
lows Eq. (|12|). Thus, for very small gap spacings (b < 0.16 
mm) and v > V we observe a growth of the interfa- 
cial width with exponent ~ 0.5 for both global and lo- 
cal scales. Similarly, in the experiments at Q = with 
b = 0.16 mm the same exponent ~ 0.65 is measured for 
both global and local scales. The power spectrum of this 
case also demonstrates the lack of dynamic scaling, as 
can be seen in the top plot of Fig. JT^(t>) . When the mod- 
ulation of the gap spacing is regular, the transition to 
saturation is sharp due to the fact that all tracks have 
the same width, and thus all the oil fingers reach the 
nominal velocity at the same time. For a random pat- 
tern the transition to saturation is smooth, and does not 
take place strictly until the widest copper tracks (those 
that have the highest value of A l + ) reach the nominal ve- 
locity. The plot for b — 0.16 mm of Fig. @ illustrates this 
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FIG. 16: Experiments with a regular modulation of the gap 
spacing, using tracks 6 mm wide. The experimental param- 
eters are v = 2V = 0.08 mm/s and b = 0.36 mm. (a) Tem- 
poral evolution of the interfaces, (b) W(l,t) plots for global 
(I = L) and local (I = L/256) scales, (c) Power spectrum at 
saturation (t = 300 s). The vertical line indicates the k cor- 
responding to the spatial periodicity of the regular pattern. 



behavior. 



VI. DISCUSSION AND CONCLUSIONS 

Our experimental results show that two ingredients are 
necessary to observe anomalous scaling in a disordered 
Hele-Shaw cell: i) strong destabilizing capillary forces, 
sufficiently persistent in the direction of growth, and 
ii) coupling in the motion between neighboring tracks. 



These two ingredients give rise to sharp slopes at the 
track edges together with a growing correlation length 
along the interface. This scenario is represented in the 
phase diagram of Fig. [[t], which shows the different re- 
gions of scaling explored experimentally. Each point of 
the diagram is an average over 6 experiments (3 dis- 
order configurations and 2 runs per disorder configura- 
tion). The dark grey area represents the region where 
an anomalous scaling has been identified. We have in- 
dicated the value of 29/ 'z, which is the value needed to 
collapse vertically the spectral curves and, therefore, in- 
dicates the degree of anomaly. For 0.15 < b < 0.60 
mm, and for driving velocities below the capil- 

lary forces are strong enough to allow an initial acceler- 
ation and subsequent deceleration of the oil fingers over 
copper tracks, but sufficiently weak to allow coupling be- 
tween neighboring tracks. For v > vm, the viscous forces 
become dominant, and the anomalous scaling is unob- 
servable. For b > 0.60 mm, the capillary forces are not 
sufficiently strong and we have ordinary scaling for any 
velocity. Finally, for b < 0.15 mm (light grey area in Fig. 
p"7| ), the capillary forces are dominant, and the oil fingers 
over copper tracks and fiberglass channels have indepen- 
dent motion. In this situation a dynamic scaling scenario 
is not adequate because there is no growing correlation 
length along the interface. Although W(l,t) grows as a 
power law of time, this is a consequence of the indepen- 
dent behavior on each track given by Eq. ([y]). Although 
it is not possible to determine the scaling exponents, still 
is possible to perform a vertical collapse of the spectral 
curves because only the difference between a and a; oc 
(i.e. 9) is needed. This allows to give the values of 29/ z 
for b = 0.16 mm. 

In conclusion, we have reported experiments of 
anomalous kinetic roughening in Hele-Shaw flows with 
quenched disorder. We have studied the displacement 
of an oil-air interface for different strength of the cap- 
illary forces, flow velocity, and disorder configurations. 
We have observed that the scaling of the interfacial fluc- 
tuations follows the intrinsic anomalous scaling ansatz 
in conditions of sufficiently strong capillary forces, low 
velocities, and a persistent disorder in the growth direc- 
tion. A random distribution of tracks has been used as 
representative of a persistent disorder. A different scaling 
for the global and local interfacial fluctuations has been 
measured, with exponents (3 ~ 0.50, (3* ~ 0.25, a ~ 1.0, 
ai oc ~ 0.5, and z ~ 2. The presence of anomalous scal- 
ing has been confirmed by the existence of multiscaling. 
When capillary forces become dominant or when a reg- 
ular track pattern is used, however, there is no growing 
correlation length along the interface, and therefore no 
dynamic scaling. A detailed investigation of the inter- 
face dynamics at the scale of the disorder shows that 
the anomaly is a consequence of the different velocities 
over copper tracks and fiberglass channels plus the cou- 
pled motion between neighboring tracks. The resulting 
dynamics gives interface profiles characterized by high 
slopes at the edges of the copper tracks. These slopes 
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follow an anomalous Levy distribution characteristic of 
systems which display anomalous scaling. 
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